Disconnected glass-glass transitions and swallowtail bifurcations 
in microscopic spin models with facilitated dynamics 
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It has been recently established that heterogeneous bootstrap percolation and related dynamic 
facilitation models exhibit a complex hierarchy of continuous and discontinuous transitions depend- 
ing on lattice connectivity and kinetic constraints. Here the range of the previously observed phase 
diagram topologies and higher-order singularities is extended to disconnected glass-glass transitions 
and swallowtail bifurcations (generic and degenerate). The phase diagram and the order parameter 
for two different types of spin mixtures are analytically determined and an experimental realization 
of the new predictions emerging in our approach is suggested. 
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Soft matter enjoys a rich variety of multiphase equi- 
libria due to the subtle interplay of energetic and en- 
tropic forces acting on different length scales. When one 
of the parameters controlling the system thermodynam- 
ics is suddenly changed, however, the phase formation is 
generally hindered for kinetic reasons and one observes 
amorphous states with distinct physical (as opposed to 
chemical) features. Such novel states of structural ar- 
rest and the glass-to-glass transition they can possibly 
undergo, were first predicted within schematic Mode- 
Coupling Theory (MCT) [l| , and have been subsequently 
observed in short-range attractive colloids, dense copoly- 
mer micellar solutions, and molecular dynamics simula- 
tions of model systems [2-12|. Their characterization is 
not only technologically important for material design 
but is also a theoretical challege, as there is no obvious 
way to discriminate on a macroscopic scale disordered 
patterns that are apparently featureless from a geomet- 
ric standpoint. 

Especially intriguing in this regard, is the existence 
of a peculiar glass-to-glass transition that has been re- 
cently identified in hard-sphere systems with a square- 
shoulder potential. In these systems the competition 
between the two repusive length scales, naturally intro- 
duced by the interaction potential, makes the glass-glass 
transition completely disconnected from the liquid phase 
in a certain range of the control parameters. In this pa- 
per I show that similar features generally exist in mi- 
croscopic on-lattice models with facilitated dynamics. A 
new prediction that emerges from our calculations is that 
the disconnected glass-glass transition can appear in sys- 
tems with either a discontinuous or a continuous liquid- 
glass transition. The framework naturally suggests that 
the latter possibility should be realized in fluid mixtures 
confined in a disordered porous matrix. Our theoretical 
analysis, which is complementary to MCT, confirms that 
these features are controlled primarily by the interplay 
of steric or kinetic hindrance effects on different micro- 
scopic lenght scales (due, e.g., to particles of dissimilar 
size), and are therefore generally expected in multicom- 
ponents systems whenever competing packing effects are 



important. 

Heterogeneous facilitation approach. In the past 
decades there has been a long lasting effort to identify 
valuable on-lattice models enabling a detailed analysis of 
the microscopic mechanism behind the glass transition. 
A promising approach in this direction is provided by the 
heterogeneous extension of bootstrap percolation and dy- 
namic facilitation ideas In this framework, the 
coarse-grained structure of a system is represented by 
an assembly of mesoscopic cells. Typically, to every cell 
i is assigned a binary spin variable, s, = ±1, depend- 
ing on the average local density (higher or lower). In 
the simplest case, no energetic interaction among cells 
is assumed, T-L = — hJ2i s i- The crucial assumption is 
that the temporal evolution of the system is dictated by 
a kinetic constraint: fluctuations in the cell i occur if 
and only if there is a certain number, say of nearby 
low-density cells. /$ is the local facilitation (or thresh- 
old) parameter which mimics the local cage effect and 
can take values in the range < < z, where z is the 
lattice connectivity. The facilitation distribution 7r(/) 
reflects the coexistence of different length scales in the 
system, e.g., mixtures of more or less mobile molecules, 
or mixtures of polymers with small and large gyration ra- 
dius. In facilitated spin mixtures the average strength of 
kinetic constraints can be tuned smoothly by changing 
the populations of spins with different fc, and one can 
thus explore a variety of different situations. Interesting 
results are obtained when the facilitated dynamics is co- 
operative, i.e., when /j > 2. Explicit calculations and 
detailed numerical simulations have shown that the ba- 
sic results of MCT schematic models are well reproduced 
within this framework [16-18]. Qualitatively, when the 
fraction of spins with /; = z — l,z is larger than that 
with 2 < fi < z — 2, the liquid- glass transition is con- 
tinuous (and the incipient cluster of frozen spins is frac- 
tal) while in the opposite case it is discontinuous (with 
a corresponding core having a compact structure). In 
the intermediate situation there is a crossover between 
the two transitions that can be either smooth [13] or 
abrupt (lij . In the latter case, the discontinuous transi- 
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tion extends deep inside the glass phase, thus generating 
an extra glass-glass transition [18j. As we are going to 
see in the remaining part of the paper, unusual features 
occur when a more subtle competition between clusters 
of frozen spins with different facilitations is present. 

Formalism. - Exploiting the analogy with heteroge- 
neous bootstrap percolation on a locally tree- like random 
graphs 15l-ll7| one can see that the probability B that a 
cell is, or can be brought, in the lower density state by 
only rearranging the state of the nearby z — 1 cells obeys 
a self-consistent polynomial equation Q(B) = where 



Q(B) = l-B-p® f z _ 1 (l-B), 
and we have defined the auxiliary function 



'(1 -X) n 




(1) 



(2) 



Here p is the fraction of high-density cells in thermal equi- 
librium at temperature T, p^ 1 = 1 + cxp(—h/k B T), and 
the angular brackets, (•••), represents the average over 
the probability distribution, 7r(/), of kinetic constraints. 
Quite generally, one finds that at high temperature ev- 
ery cell can always change state, B — 1, while at low 
temperature there is a fraction of spins which is unable 
to change state, B < 1, and the system is therefore a 
glass. The detailed topology of the phase diagram at low 
temperature depends on the coefficients of the polyno- 
mial Q{B) and can be rather intricate. Interestingly, the 
self-consistent equation Q(B) = has a formal structure 
quite similar to that satisfied by the long-time limit of 
the form factor in MCT. Accordingly, one can immedi- 
ately draw the conclusion that Q(B) exhibits the same 
hierarchy of bifurcations, Ai, of schematic MCT. They 
occur when the corresponding maximum root of Q has a 
degeneracy t (£ > 2) and 



dB n 



= 0, n = 0,. ..,£-!; 



(3) 



The Taylor expansion of Q near the critical surface and 
Eqs. ©, immediately implies that the scaling form of the 
order parameter near an Ai bifurcation goes like e 1 ^, 
where e is the distance from the critical surface (e.g., 
e = T—T c ). Singularities of type Ai can be further distin- 
guished in generic and degenerate depending on whether 
the order parameter, changes abruptly or smoothly 
near the transition. To denote this latter case, we shall 
use the notation A|. 

In the following, we shall focus on ternary mixtures 
with facilitation distribution 



"■(/i) = (! - l)3fi,a + (q- r)Sf ilb + rS fi 



(4) 



For such ternary mixtures, denoted here with (a, b, c), 
we shall consider two distinct situations corresponding to 



small and large facilitation variance per spin population 
o-j = ((/ 2 ) — (f) 2 )/3. For each situation, we determine 
the phase diagram and compute the fraction of perma- 
nently frozen spins, $, which represents the actual order 
parameter in this framework. It is directly related to B 
through the general relation 



$ =p *f(B) + (1 -pW,{p%Z{{B)) 



(5) 



For sake of simplicity we shall consider hereafter only 
random graphs with fixed connectivity, i.e. Bethe lat- 
tices. Similar results are expected in more general ran- 
dom graphs with variable connectivity 16, 19, 2p| |. 




q — r 

FIG. 1. A section of the phase diagram for the mixture 
(2, 3, 4) on a Bethe lattice with connectivity z = 5. Thick lines 
represent discontinuous liquid-glass and glass-glass transitions 
while thin dashed lines correspond to continuous liquid-glass 
transition. For r — 0.34 the discontinuous glass-glass transi- 
tion is disconnected from the liquid phase. 



Mixture with small a 2 - . Let us first consider the 
(2,3,4) mixture on a Bethe lattice with z = 5. For such 
a mixture, the fixed-point equation obeyed by B is: 



1 

P 



l + B 



3B 3 + 6qB 2 {l - B) +4rB 3 . (6) 



The continuous glass transition, T c (r) = — l/ln(4r — 1), 
is obtained by setting B = 1 in Eq. ©. Therefore, it 
does not depend on q and is limited to the range 1/2 > 
r > 1/4 (we do not consider here the case of negative 
temperature). Setting the first-order derivative of Eq. © 
is zero, we get 



-9B 2 + WB - 1 - YlrB 2 
12B - 18B 2 



(7) 



and thus the discontinuous transition is obtained by plot- 
ting Eqs. ((6|) and Q parametrically in terms of B. A sec- 
tion of the phase diagram illustrating the distinct topolo- 
gies of the transition lines obtained for different values of 
r (i.e., the fraction of spins with / = 4), is shown in the 
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FIG. 2. Higher-order singularities in the mixture (2, 3, 4) on 
a Bethe lattice with z = 5. The upper and lower curves 
represent the line of generic and degenerate cusp bifurcation, 
respectively. The two curves touches smoothly at a degener- 
ate AJ swallowtail bifurcation. The dashed and dotted curves 
belong to the unstable branches of the phase diagram. 



Fig. Q] For small values of r we see that there is only 
a discontinuous liquid-glass transition (corresponding to 
the fold bifurcation). The continuous transition occurs 
for r > 0.27 and crosses the discontinuous line. This 
latter, in turn, enters the frozen phase thus producing 
an extra glass-glass transition. Upon further increasing 
r something more interesting happens: we observe that 
the glass-to-glass transition departs from the continuous 
liquid-glass transition and becomes completely discon- 
nected from the liquid phase. This departure generates 
an extra endpoint singularity corresponding to a generic 
cusp bifurcation, A 3 . The range of q values over which 
the transition is disconnected widens until the the glass- 
glass transition eventually disappears. 

The endpoints of the glass-glass transition lines de- 
fine a set of generic cusp singularities, A3, whereas the 
crossing points between the continuous and discontinu- 
ous liquid-glass transitions correspond to degenerate cusp 
singularities, A3. They are respectively given by 
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ln(15 - ISq), 
T(q) ~ { -ln(2g-l), 



A 3 ; 
At 



(8) 



Either curves possess an unstable branch and are repre- 
sented in the (T, q) plane in Fig. [2] One can easily check 
that they coaelesce smoothly in a degenerate swallow- 
tail singularity, A|, which is exactly located at q = 2/3, 
r = 1/3, and T = l/ln3. Notice that these values of q 
and r corresponds to a mixture with an exactly balanced 
composition of each components (1— q = q—r = r = 1/3). 
For completeness we also show, in Fig.[3j the variation of 
the order parameter with for a value of r in the range in 
which we observe the disconnected glass-glass transition. 

The above findings have not been previously reported 
and we expect they should be observed in other systems, 



FIG. 3. The order parameter, representing the fraction of 
frozen spins vs temperature, T, at r = 0.32 and for several 
values of q in a mixture (2, 3, 4) with 2 = 5. 



e.g., in fluid mixtures confined in porous media. In fact, 
building on the observations of Refs. (21, 22 1, spin-glass 
models with p + q multispin interaction terms when sup- 
plemented with an extra random field should reproduce 
the structure of MCT for a fluid binary mixture in a ran- 
dom environment. However, the kinetic analogue of the 
three Hamiltonian interaction terms correspond just to 
the facilitated terms entering the model studied above. 



Mixture with large er^ . - Next we consider, for com- 
parison, the (2,5,7) mixture with probability distribu- 
tion Eq. (0| on a Bethe lattice with z = 12. Evidently, 
since there is no spin population with = z — 1, z, 
the order parameter cannot vanish continuously and so 
no degenerate singularity is expected. A section of the 
phase diagram obtained in this case is reported in the 
Fig. [4] for several values of r illustrating the different 
topologies of the transition lines we obtain. For r < 0.4 
we find that there is only a discontinuous liquid-glass 
transition, i.e., a fold bifurcation, A2 singularity. For 
r > 0.4 it appears an additional cusp bifurcation, A 3 , 
corresponding to the end point of a glass-glass transi- 
tion line. The crossover between these two behaviors oc- 
curs via a swallowtail bifurcation, A4 , which is located at 
q ~ 0.773, T ~ 0.296, r ~ 0.40369. Upon increasing fur- 
ther r the glass-glass transition becomes eventually dis- 
connected from the liquid phase, see the curve r = 0.68 in 
Fig.|H Therefore, also in this case, spin populations with 
different values of fi, with 2 < fa < z — 2, compete with 
each other to produce a disconnected glass-glass transi- 
tion. As we can observe in Fig. 2] this occurs when the 
fraction 1 — q of spins with / = 2 become smaller than 
that with f — 7 (i.e., r) and, correspondingly, the frac- 
tion of spins with intermediate value of /, / = 5, becomes 
pretty small, q « r. The variation of the order parame- 
ter $ with the temperature, T, is shown in Fig. [5] for r 
in the range in which we observe the disconnected glass- 
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FIG. 4. A section of the phase diagram for the mixture / = 
(2, 5, 7) on a Bethe lattice with connectivity z = 12. The lines 
represent discontinuous transitions. In this figure, the glass- 
glass transition occurs for r = 052, 0.6 and 0.68. In the latter 
case it is disconnected from the liquid phase. Higher-order 
singularities with cusp and swallowtail structures are denoted 
with A3 and A4 respectively. 
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FIG. 5. The order parameter, <3>, representing the fraction 
of frozen spins vs temperature, T, at r = 0.7 and for several 
values of q, in the (2, 5, 7) mixture with z — 12. The double 
jump of "I" for q = 0.74 and 0.76 corresponds to a liquid-glass 
transition followed by a glass-glass transition upon lowering 
the temperature. 



glass transition and several values of q. As expected we 
find that exhibits a single or a double jumps depending 
on whether the system crosses one or two transition lines 
upon lowering the temperature. 

The results we find here are qualitatively similar to 
those obtained in hard-sphere systems with the square- 
shoulder potential [ll|. Nevertheless, in the latter case, 
there is an interesting extra feature which is not re- 
produced in our approach: It is the counter-intuitive 
melting-by-cooling often associated with higher-order 
singularities. It would be interesting to ascertain the in- 
gredients that are actually required to observe this reen- 
trancy effect in the present context, whether they are 
entropic or energetic in origin. 

It is also interesting to notice that the phase diagram 
we have derived above shares some similarities with that 
of the spherical p + q spin-glass models (see Fig. 1 in 
Ref. The disconnected glass-glass transition, how- 

ever, is not observed in this spin-glass model. Since a 
p— spin interaction term plays the same role of a facili- 
tated dynamics with z — 2 > /, > 2, our approach sug- 
gests that an additional multispin interaction term in the 
p + q Hamiltonian is required to observe a disconnected 
glass-glass transition. Nevertheless, the analysis of this 
latter case could be rather awkward as it involves delicate 
aspects of replica symmetry breaking calculations 

Conclusions.- To summarise, I have shown that dis- 
connected glass-glass transitions take place in systems 
with either a continuous or a discontinuous liquid-glass 
transition and that the related swallowtail bifurcation 
can be degenerate or generic, respectively. Such fea- 
tures are controlled primarily by the competition between 



packing effects on different microscopic length scales, as 
opposed to system-specific details of the molecular inter- 
actions, and are therefore expected in a range of soft mat- 
ter systems including fluid mixtures confined in porous 
media, and colloidal and polymer gels. The fact that, 
starting from very different premises, we reach conclu- 
sions quite similar to those obtained in more realistic 
model systems is neither obvious nor coincidental, and 
should lend further support to the universality of both 
MCT and the present framework. 

Although the dynamical behavior of facilitated systems 
can be numerically studied with a relatively modest effort 
(using continuous-time algorithms) it would be certainly 
interesting, from a theoretical perspective, to rationalise 
the anomalous logarithmic relaxation near higher-order 
singularities through the dynamics of minimal size rear- 
rangements 25|. Similarly, one should be able to iden- 



tify the conditions under which the parameter exponent, 
typically denoted with A in MCT, can be derived in the 
present starting from the "facilitated" analogue of the 
MCT functional That would provide an alternative 
route to the determination of A by-passing the computa- 
tion of multispin correlation functions recently suggested 
in Ref. |26| . Work in these directions is in progress and 
we hope to report on such issues in the near future. 
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